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We consider the construction of generic spherically symmetric thin-shell traversable wormhole 
, spacetimes in standard general relativity. By using the cut-and-paste procedure, we comprehensively 

• analyze the stability of arbitrary spherically symmetric thin-shell wormholes to linearized spherically 

symmetric perturbations around static solutions. While a number of special cases have previously 
been dealt with in scattered parts of the literature, herein we take considerable effort to make the 
analysis as general and unified as practicable. We demonstrate in full generality that stability of 
the wormhole is equivalent to choosing suitable properties for the exotic material residing on the 
wormhole throat. 
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PsJ ' Traversable wormholes are hypothetical tunnels in spacetime, through which observers may freely travel [1, 2]. 
. These geometries are supported by "exotic matter" , involving a stress-energy tensor violating the null energy condition 
■ (NEC). That is, there exists at least one null vector fc'' such that T^^k'^'-k'^ < on or in the immediate vicinity of the 
T-H . wormhole throat (see, in particular, [3-5]). In fact, wormhole geometries violate all the standard pointwise energy 
I ■ conditions, and all the standard averaged energy conditions [3] . Although (most) classical forms of matter are believed 
^ ] to obey (most of) the standard energy conditions [6] , it is a well-known fact that they are violated by certain quantum 
• ^ ■ effects, amongst which we may refer to the Casimir effect, Hawking evaporation, and vacuum polarization. (See [3], 
I and more recently [7], for a review. For some technical details see [8].) It is interesting to note that the known 
violations of the pointwise energy conditions led researchers to consider the possibility of averaging of the energy 
conditions over timelike or null geodesies [9]. For instance, the averaged weak energy condition (AWEC) states that 
the integral of the energy density measured by a geodesic observer is non-negative. (That is, / T^^yU^^U" dr > 0, 
where r is the observer's proper time.) Thus, the averaged energy conditions are weaker than the pointwise energy 
conditions, they permit localized violations of the energy conditions, as long as they hold when suitably averaged 
along a null or timelike geodesic [9] . 

As the theme of exotic matter is a problematic issue, it is useful to minimize its usage. In fact, it is important to 
emphasize that the theorems which guarantee the energy condition violation are remarkably silent when it comes to 
making quantitative statements regarding the "total amount" of energy condition violating matter in the spacetime. 
In this context, a suitable measure for quantifying this notion was developed in [10], see also [11], where it was 
shown that wormhole geometries may in principle be supported by arbitrarily small quantities of exotic matter (an 
interesting application of the quantification of the total amount of energy condition violating matter in warp drive 
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spacetimes was considered in [12]). In the context of minimizing the usage of exotic matter, it was also found that 
for specific models of stationary and axially symmetric traversable wormholes the exotic matter is confined to certain 
regions around the wormhole throat, so that certain classes of geodesies traversing the wormhole need not encounter 
any energy condition violating matter [13]. For dynamic wormholes the null energy condition, more precisely the 
averaged null energy condition, can be avoided in certain regions [14-16]. Evolving wormhole geometries were also 
found which exhibit "flashes" of weak energy condition (WEC) violation, where the matter threading the wormhole 
violates the energy conditions for small intervals of time [17]. In the context of nonlinear electrodynamics, it was 
found that certain dynamic wormhole solutions obey (suitably defined versions of) the WEC [18]. 

It is interesting to note that in modified theories of gravity, more specifically in f{R) gravity, the matter threading 
the wormhole throat can be forced to obey (suitably defined versions of) all the energy conditions, and it is the 
higher-order curvature terms that are responsible for supporting these wormhole geometries [19]. (Related issues arise 
when scalar fields are conformally coupled to gravity [20, 21].) In brancworlds, it was found that it is a combination of 
the local high-energy bulk effects, and the nonlocal corrections from the Weyl curvature in the bulk, that may induce 
an effective NEC violating signature on the brane, while the "physical" / "observable" stress-energy tensor confined to 
the brane, threading the wormhole throat, nevertheless satisfies the energy conditions [22]. 

An interesting and efficient manner to minimize the violation of the null energy condition, extensively analyzed in 
the literature, is to construct thin-shell wormholes using the thin-shell formalism [3, 23, 24] and the cut-and-paste 
procedure as described in [3, 25-28]. Motivated in minimizing the usage of exotic matter, the thin-shell construction 
was generalized to nonspherically symmetric cases [3, 25], and in particular, it was found that a traveler may traverse 
through such a wormhole without encountering regions of exotic matter. In the context of a (limited) stability analysis, 
in [26] , two Schwarzschild spacetimes were surgically grafted together in such a way that no event horizon is permitted 
to form. This surgery concentrates a nonzero stress energy on the boundary layer between the two asymptotically fiat 
regions and a dynamical stability analysis (with respect to spherically symmetric perturbations) was explored. In the 
latter stability analysis, constraints were found on the equation of state of the exotic matter that comprises the throat 
of the wormhole. Indeed, the stability of the latter thin-shell wormholes was considered for certain specially chosen 
equations of state [3, 26], where the analysis addressed the issue of stability in the sense of proving bounded motion 
for the wormhole throat. This dynamical analysis was generalized to the stability of spherically symmetric thin-shell 
wormholes by considering linearized radial perturbations around some assumed static solution of the Einstein field 
equations, without the need to specify an equation of state [28]. This linearized stability analysis around a static 
solution was soon generalized to the presence of charge [29] , and of a cosmological constant [30] , and was subsequently 
extended to a plethora of individual scenarios [31, 32], some of them rather ad hoc. 

The key point of the present paper is to develop an extremely general, flexible, and robust framework that can 
quickly be adapted to general spherically symmetric traversable wormholes in 3-1-1 dimensions. We shall consider 
standard general relativity, with traversable wormholes that are spherically symmetric, with all of the exotic material 
confined to a thin shell. The hulk spacetimes on either side of the wormhole throat will be spherically symmetric 
and static but otherwise arbitrary (so the formalism is simultaneously capable of dealing with wormholes embedded 
in Schwarzschild, Reissner-Nordstrom, Kottler, or de Sitter spacetimes, or even "stringy" black hole spacetimes). 
The thin shell (wormhole throat), while constrained by spherical symmetry, will otherwise be permitted to move 
freely in the bulk spacetimes, permitting a fully dynamic analysis. This will then allow us to perform a general 
stability analysis against spherically symmetric perturbations, where wormhole stability is related to the properties 
of the exotic matter residing on the wormhole throat. We particularly emphasize that our analysis can deal with 
geometrically imposed "external forces", (to be more fully explained below), a feature that has so far been missing 
from the published literature. Additionally we emphasize the derivation of rather explicit and very general rules 
relating the internal structure of the wormhole throat to a "potential" that drives the motion of the throat. 

This paper is organized in the following manner: In Sec. II we outline in detail the general formalism of generic 
dynamic spherically symmetric thin-shell wormholes, and provide a novel approach to the linearized stability analysis 
around a static solution. In Sec. Ill, we provide specific examples by applying the generic linearized stability formalism 
outlined in the previous section. In Sec. IV, we draw some general conclusions. 

II. FORMALISM 

Let us first perform a general theoretical analysis; subsequently we shall look at a number of specific examples. 
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A. Bulk spacetimes 

To set the stage, consider two distinct spacetime manifolds, and A^_, with metrics given by g'^^{x'^) and 
.9/7i/(2^-)i i'^ terms of independently defined coordinate systems x'^ and x^. A single manifold M is obtained by gluing 
together the two distinct manifolds, and Ai-, i.e., J\4 — U Ai-, at their boundaries. The latter are given 
by and E_, respectively, with the natural identification of the boundaries E = S+ = E_. This construction is 
depicted in the embedding diagram in Fig. 1, and is further explored in the sections given below. 




FIG. 1: The figure depicts an embedding diagram of a traversable thin-shell wormhole. A single manifold Ml is obtained by 
gluing together two distinct spacetime manifolds, M+ and M-, at their boundaries, i.e., M — U M-, with the natural 
identification of the boundaries E = E+ = E_. 



More specifically, throughout this work, we consider two generic static spherically symmetric spacetimes given by 
the following line elements: 

-1 





dt\ + 











dr\ 



(1) 



on 7W±, respectively. Using the Einstein field equation, G^^ — SirTf^^, (with c = G = 1), the (orthonormal) stress- 
energy tensor components in the bulk are given by 

= z^b', (2) 



Pr{r) 
Pt{r) 



87rr2 
1 



Sirr 
1 



b', 

-[2$' (6 



:[{-b + 3rb' - 2r)$' + 2r{b - r)($') + 2r{b - r)$" + b"r] 



(3) 
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where the prime denotes a derivative with respect to the radial coordinate. Here p{r) is the energy density, Pr{r) is 
the radial pressure, and pt (r) is the lateral pressure measured in the orthogonal direction to the radial direction. The 
± subscripts were (temporarily) dropped so as not to overload the notation. 

The energy conditions will play an important role in the analysis that follows, so we will at this stage define the 
NEC. The latter is satisfied if r^,^ k^^ k'^ > 0, where T^i, is the stress-energy tensor and fc^ any null vector. Along the 
radial direction, with fc'' — (1, ±1,0,0) in the orthonormal frame, where Tp,c, = dia.g[p{r),pr{r),pt{r),pt{r)], we then 
have the particularly simple condition 



p{r) +Pr{r) 



(r-6)$ 
47rr2 



> 0. 



(5) 
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Note that in any region where the t coordinate is timehke (requiring r > 6(r)) the radial NEC reduces to $'(r) > 0. 
The NEC in the transverse direction, p + pt > 0, does not have any direct simple interpretation in terms of the metric 
components. 

We emphasize that the results outlined in this work are also valid for the case of an intra- universe wormhole, i.e., 
with a single manifold with a wormhole connecting distant regions. This is indeed true as long as the bulk geometries 
are both asymptotically flat, then they can be viewed as widely separated parts of the same asymptotically flat 
spacetime (to an arbitrarily good approximation that gets better as the two wormhole mouths get further and further 
separated). For instance, see Fig. 2 for a spacetime diagram depicting a thin-shell wormhole in an asymptotically flat 
spacetime, represented by two manifolds joined at a radial coordinate a. 



future timelike infinity 




past timelike intinity 



FIG. 2: The spacetime diagram for a thin-shell wormhole in an asymptotically flat spacetime, represented by two manifolds 
M.+ and A4-, joined at a radial coordinate a, at the junction surface S. 



B. Extrinsic curvature 

The manifolds are bounded by hypersurfaces and respectively, with induced metrics g^j and g^j. The 
hypersurfaces are isometric, i.e., g^j{£,) — 5ij(C) — dijiOj i'^ terms of the intrinsic coordinates, invariant under the 
isometry. As mentioned above, single manifold Ai is obtained by gluing together and Ai- at their boundaries, 
i.e., Ai — Ai+ U A4^, with the natural identification of the boundaries S — S+ — S_. The three holonomic basis 
vectors e^j) = d/d^^ tangent to S have the following components e|'j^|± = 9x^/9^*, which provide the induced metric 

on the junction surface by the following scalar product gij = ej,;) • e^j) = .9/^iye^--)eQ-^ |±. The intrinsic metric to E is 
thus provided by 

dsl ^ -dr^ +a(Tf {d0^ +ain'^9d(t>'^). (6) 

Thus, for the static and spherically symmetric spacetime considered in this work, the single manifold, A4, is 
obtained by gluing M+ and M- at S, i.e., at f{r,T) = r — a{T) = 0. The position of the junction surface is given 
by x'^{t, 9, (j)) = {t{T) , a{T) , 9 , (fi), and the respective 4- velocities (as measured in the static coordinate systems on the 
two sides of the junction) are 





where the overdot denotes a derivative with respect to r, which is the proper time of an observer comoving with the 
junction surface. 
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Wc shall consider a timelike junction surface S, defined by the parametric equation of the form f{x^{^^)) = 0. The 
unit normal 4— vector, n^, to S is defined as 



± 



-1/2 



with = +1 and n^eJ^^s^ = 0. The Israel formalism requires that the normals point from A4- to A4+ [24]. Thus, 

the unit normals to the junction surface, determined by Eq. (8), are given by 



,,™U^,0,0 . (9) 



^ _ b±{a) V a 



Note (in view of the spherical symmetry) that the above expressions can also be deduced from the contractions 
J/'^n^ = and n^n^ = +1. The extrinsic curvature, or the second fundamental form, is defined as Kij = ni^.^^e'^.-^e'^jy 

Differentiating n^e^.^ = with respect to we have 



d'^x'' dx^^dx" 

= -n^,^l^^, (10) 



so that in general the extrinsic curvature is given by 

Note that for the case of a thin shell Kij is not continuous across E, so that for notational convenience, the discontinuity 
in the second fundamental form is defined as kh = Kf'- — K~-. 

Using equation (11), the non-trivial components of the extrinsic curvature can easily be computed to be 



= +'^'^ia)Jl-^-^ + a^}, (13) 

where the prime now denotes a derivative with respect to the coordinate a. 
Several comments are in order: 

• Note that is independent of the quantities a circumstance that will have implications later in this 
article. This is most easily verified by noting that in terms of the normal distance £ to the shell S the extrinsic 
curvature can be written as Kij — \digij — ^n'^d^gij = ^n^drgij, where the last step relies on the fact that the 
bulk spacetimes are static. Then since goo = r^, differentiating and setting r — >■ a we have Kgg = arf . Finally 

= ^, (14) 

which is a particularly simple formula in terms of the radial component of the normal vector, and which easily 

lets us verify (12). 

• For Ktt there is an argument (easily extendable to the present context) in Ref. [3] (see especially pages 181-183) 
to the eflfect that 

K"^^ = ± (magnitude of the physical 4-acceleration of the throat). (15) 
This gives a clear physical interpretation to K'^^ and rapidly allows one to verify (13). 

• There is also an important differential relationship between these extrinsic curvature components 

^{ae''^ K%^} = e^^ K%^a. (16) 
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The most direct way to verify this is to simply differentiate, using (12) and (13) above. Geometricahy, the 
existence of these relations between the extrinsic curvature components is ultimately due to the fact that the 
bulk spacetimes have been chosen to be static. By noting that 

we can also write this differential relation as 

{a e*± K^g^} = e*± K\^. (18) 

• We emphasize that including the possibility that $± (a) ^ is already a significant generalization of the extant 
literature. 



C. Lanczos equations: Surface stress-energy 

The Lanczos equations follow from the Einstein equations applied to the hypersurface joining the bulk spacetimes, 
and are given by 



(19) 



Here S'^j is the surface stress-energy tensor on S. In particular, because of spherical symmetry considerable simplifica- 
tions occur, namely n^j = diag (k"^^, k^^, k^q). The surface stress-energy tensor may be written in terms of the surface 
energy density, cr, and the surface pressure, V, as S^^ = dia,g{—a,V,V). The Lanczos equations then reduce to 



(20) 
(21) 



Taking into account the computed extrinsic curvatures, Eqs. (12)-(13), we see that Eqs. (20)-(21) provide us with 
the following expressions for the surface stresses: 



1 

47ra 



1 - 



h+{a) 



(22) 



1 

Stto 



1 - 



6+ (a) 



+ 



l + a2 



h- {a)-\-ah'_ (a) 
2a 



1 



b-{a) 



^2 



6_(a) 



+ b? a$'_(a) 



(23) 



The surface stress-energy tensor on the junction surface S is depicted in the embedding spacetime diagram in Fig. 3, 
in terms of the surface energy density, a, and the surface pressure, V. 

Note that the surface mass of the thin shell is given by rUg = ^-KO^a, a quantity which will be extensively used below. 
Note further that the surface energy density cr is always negative, (which is where the energy condition violations show 
up in this thin-shell context), and furthermore that the surface energy density a is independent of the two quantities 



D. Conservation identity 

The first contracted Gauss-Codazzi equation, sometimes called simply the Gauss equation, or in general relativity 
more often referred to as the "Hamiltonian constraint" , is 

G^. n-^ = i [K^ - Ki,K'^ - ^i?) . (24) 
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FIG. 3: The figure depicts an embedding diagram of a traversable tliin-siiell wormliole. The single manifold A4 is obtained by 
gluing together M+ and A^_, at junction surface E. The surface stress-energy tensor on E is given in terms of the surface 
energy density, a, and the surface pressure, V. See the text for details. 



Together with the Einstein equations this provides the evolution identity 



(25) 



The convention [X\^ = — and X = ^{X^\y: + -'^"Is) is used. The second contracted Gauss-Codazzi 

equation, sometimes called simply the Codazzi or the Codazzi-Mainardi equation, or in general relativity more often 
referred to as the "ADM constraint" or "momentum constraint" , is 



Together with the Lanczos equations this provides the conservation identity 



(26) 



^ 3\i - 



(27) 



When interpreting the conservation identity Eq. (27), consider first the momentum flux defined by 

+ 



1 + 



±{Tii + Tff) 



1 



(28) 



where Tff and Tff are the bulk stress-energy tensor components given in an orthonormal basis. This flux term 
corresponds to the net discontinuity in the (bulk) momentum flux = T^j, which impinges on the shell. (This 
flux term is identically zero in all the currently extant literature.) Applying the (bulk) Einstein equations we see 



It is useful to deflne the quantity 



a 
Ana 



$;(«)./ 1 



h + 5>_(a)y 1 h 



1 



(29) 



(30) 
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and to let A = Atto^ be the surface area of the thin shell. Then in the general case, the conservation identity provides 
the following relationship 



or equivalently 



da , ^^ 1 dA „ . 



d(aA) ^dA „ , . , „, 

^+P-=S^a, (32) 



The first term represents the variation of the internal energy of the shell, the second term is the work done by the 
shell's internal force, and the third term represents the work done by the external forces. Once could also brute force 
verify this equation by explicitly differentiating (22) using (23) and the relations (16). If we assume that the equations 
of motion can be integrated to determine the surface energy density as a function of radius a, that is, assuming the 
existence of a suitable function (7(a), then the conservation equation can be written as 

a' = ~{a + r) + E, (33) 

where a' = da /da. We shall carefully analyze the intcgrability conditions for a {a) in the next sub-section. For 
now, note that the flux term (external force term) S is automatically zero whenever $± = 0; this is actually a 
quite common occurrence, for instance in either Schwarzschild or Reissner-Nordstrom geometries, or more generally 
whenever p + pr = 0, so it is very easy for one to be mislead by those special cases. In particular, in situations of 
vanishing flux S = one obtains the so-called "transparency condition". [G^,y ?7^n'^]^ = 0, see [32]. The conservation 
identity, Eq. (27), then reduces to the simple relationship a = —2{a + V)a/a. But in general the "transparency 
condition" does not hold, and one needs the full version of the conservation equation as given in Eq. (31). 

Physically, it is particularly important to realize that merely specifying the two bulk geometries [g±]a6 is not enough 
to fully determine the motion of the throat — some sort of assumption must be made regarding the internal behaviour 
of the physical material concentrated on the throat itself, and we now turn to investigating this critically important 
point. 



E. Integrability of the surface energy density 

When does it make sense to assert the existence of a function cr(a)? Let us start with the situation in the absence 
of external forces (we will rapidly generalize this) where the conservation equation, 

a = -2{a + V)ala (34) 

can easily be rearranged to 

" =-2^. (35) 



a + V a 

Assuming a barotropic equation of state 'P{a) for the matter on the wormhole throat, this can be integrated to yield 

d(T ^ da „ , / , X ,„„s 

-2 1 — = -21n(a/ao). (36) 



cr + V{a) a 

This implies that a can be given as some function a(cr) of a, and by the inverse function theorem implies over suitable 
domains the existence of a function a{a). Now this barotropic equation of state is a rather strong assumption, albeit 
one that is very often implicitly made when dealing with thin-shell wormholes (or thin-shell gravastars [33-35], or 
other thin-shell objects). As a first generalization, consider what happens if the surface pressure is generalized to be 
of the form 'P{a, a), which is not barotropic. Then the conservation equation can be rearranged to be 

^,^_2[a + na,a)]_ ^3^^ 



This is a first-order (albeit nonlinear and non-autonomous) ordinary differential equation, which at least locally will 

have solutions a (a). There is no particular reason to be concerned about the question of global solutions to this ODE, 
since in applications one is most typically dealing with linearization around a static solution. 
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If we now switch on external forces, that is <J>^ 7^ 0, then one way of guaranteeing integrabihty would be to demand 
that the external forces are of the form R{a,a), since then the conservation equation would read 

a' = -'±±^+Eia,a), (38) 

which is again a first-order albeit nonlinear and non-autonomous ordinary differential equation. But how general is 
this S = S(a, a) assumption? There are at least two situations where this definitely holds: 

• If = ^-(a) = ^{d), then S = — $'(a) a, which is explicitly of the required form. 

• If 6+(a) = b-{a) = 6(a), then S = if^YC*^) + ^'-{d)] cr, which is again explicitly of the required form. 

But in general we will need a more complicated set of assiimptions to assure intcgrability, and the consequent existence 
of some function a{a). A model that is always sufficient (not necessary) to guarantee integrability is to view the 
exotic material on the throat as a two-fluid system, characterized by a± and 'P±, with two (possibly independent) 
equations of state V±{(j±). Specifically, take 

C7± = -^{K±)%, (39) 
■P± = ^{{K±r, + {K±)%} . (40) 



In view of the differential identities 



dr 

each of these two fluids is independently subject to 



^ {a e'^^ K%^} = e"^^ a, (41) 



^{e*±a±} = -^{a±+P±}d. (42) 

which is equivalent to 

{e*±a±}' = -^{a±+P±}. (43) 

With two equations of state 'P±{a±) these are two nonlinear first-order ordinary differential equations for a±. These 
equations are integrable, implicitly defining functions cr±{a), at least locally. Once this is done we define 

a{a) = a+ (a) + (t_ (a) , (44) 

and 

ms{a) = 47rcr(a) a^. (45) 

While the argument is more complicated than one might have expected, the end result is easy to interpret: We can 

simply choose cr(a), or equivalently ms{a), as an arbitrarily specifiable function that encodes the (otherwise unknown) 
physics of the specific form of exotic matter residing on the wormhole throat. 

F. Equation of motion 

To qualitatively analyze the stability of the wormhole, assuming integrability of the surface energy density, (that 
is, the existence of a function a'(a)), it is useful to rearrange Eq. (22) into the form 



1 - + = -^1 - ^ + - 4.a a{a) . (46) 

Squaring this equation, rearranging terms to isolate the square root, and finally by squaring once again, we deduce 
the thin-shell equation of motion given by 

^a^ + V{a)^0, (47) 
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where the potential V{a) is given by 





'ms{ay 


1 


r A(a) 1 




2a 




rrisia) 





(48) 



Here ms(a) = Amd^ a{a) is the mass of the thin shell. The quantities b{a) and A(a) are defined, for simplicity, as 



6(a) 
A(a) 



b+{a)+b-{a) 
2 

6+ (a) — b-{a) 



(49) 
(50) 



respectively. This gives the potential V{a) as a function of the surface mass ms{a). By differentiating with respect 
to a, (using (17)), we see that the equation of motion implies 



-V'{a). 



(51) 



It is sometimes useful to reverse the logic flow and determine the surface mass as a function of the potential. Following 
the techniques used in [34], suitably modified for the present wormhole context, a brief calculation yields 



miia) = 2a^ <( 1 - M _ 2V{a) + 



/l_?M_4na) 



2V{a) + 



b{a) 



a2 J ' 



(52) 



where one is forced to take the positive root to guarantee reality of the surface mass. Noting that the radical factorizes 
we see 



ml{a) = 



1 _ M _ 2V{a) 



l_M^_2y(a)Jl-^ 



2Y{a) 



and in fact 



ms{a) 



6+ (a) 



2V{a) 



b-ia) 



2V{a) 



(53) 



(54) 



with the negative root now being necessary for compatibility with the Lanczos equations. Note the logic here — 
assuming intcgrabiUty of the surface energy density, if we want a specific V{a) this tells us how much surface mass 
we need to put on the wormhole throat (as a function of a), which is implicitly making demands on the equation of 
state of the exotic matter residing on the wormhole throat. In a completely analogous manner, the assumption of 
integrability of a{a) implies that after imposing' the (Yjuation of motion for thc^ shc41 one has 



cr(a) = --^ 
4-n-a 



1 _ _ 2^(^) 



1 - ^ - 2V{a) 



(55) 



while 



1 - 2Via) - aVja) - 



+ \ 1 



b+ia) 



21/ (a) a$+(o) 



+ 



1 - 2V{a) - aV'{a) 



b— {a)+ab'_ (a) 



2a 



b-{a) 



2V{a) 



6_(a) 



2V{a) a$'_(a) 



and 



S(a) 



47ra 



$V(a)Wl 



b+{a) 



2V{a) + ^'_{a)\ 1 



6_(a) 



2V{a) 



(56) 



(57) 



The three quantities {(j{a), ^{a), ^{a)} (or cquivalently {ms{a), V{a), S(a)}) are related by the differential conser- 
vation law, so at most two of them are functionally independent. 
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G. Linearized equation of motion 

Consider a linearization around an assumed static solution (at ao) to the equation of motion ^d^ + V{a) = 0, and 
so also a solution of d = —V'{a). Generally a Taylor expansion of V{a) around ao to second order yields 

V{a) = V{ao) + V'{ao){a - ao) + ^V^"(ao)(a - ao)" + 0[(a - ao)'] . (58) 

But since we are expanding around a static solution, do = do = 0, we automatically have V(ao) = V'{ao) = 0, so it is 
sufHcient to consider 



Via) = \v"{ao){a - ao)" + 0[{a - ao)'] . 



(59) 



The assumed static solution at ao is stable if and only if V{a) has a local minimum at ao, which requires y"(ao) > 0. 
This will be our primary criterion for wormhole stability, though it will be useful to rephrase it in terms of more basic 
quantities. 

For instance, it is extremely useful to express m'g{a) and m"{a) by the following expressions: 



m'^{a) = + 



m,(a) a J {b+{a)/ay + 2V'{a) 



+ 



{b-{a)/ay + 2V'{a) 
2 \^/l-b+{a)/a-2V{a) ' ^yl - b-{a)/a - 2V{a) 



+ 



(60) 



and 



m'Ua) = 



{b+{a)/ay + 2V'ia) 



{b^{a)/ay + 2V'{a) 



,yi-h+{a)/a-2V{a) ^1 - b^{a)/a - 2V{a) 



+- 



[(6+(a)/a)' + 2y'(a)]^ 



4 \ [1 - b+{a)/a - 2V{a)]^/^ ^ [1 - 6_(a)/a - 2V{a)]^/^ / 
(6+(a)/a)" + 2y"(a) ^ (6_(a)/a)" + 2F"(a) 



[ib.ia)/ay + 2V'{a)]^ \ 



2 [^^l-b+{a)/a-2V{a) ^/l - b-{a)/a - 2V{a) j ' 



(61) 



Doing so allows us to easily study linearized stability, and to develop a simple inequality on m"(ao) by using the 

constraint V"(ao) > 0. Similar formulae hold for <j'{a), <T"{a), for Via), 'P"{a), and for S'(a), S"(a). In view of the 
redundancies coming from the relations ms{a) = 47r (T(a)a" and the difFcrcntial conservation law, the only interesting 
quantities are S'(a), S"(a). 

For practical calculations, it is extremely useful to consider the dimensionless quantity 



ms{a) 



= 47rcr(a)a = — 



b+{a) 



2V{a) + dl 



b-{a) 



- 2V{a) 



(62) 



and then to express [ms{a)/a]' and [ms(a)/o]" by the following expressions: 



and 



[ms{a)/ay = +- 



[ms{a)/a\" = +- 



1 

+ 2 



It is similarly useful to consider 



47rS(a)a = 



{b+{a)/ay + 2V'{a) 



{b.{a)/ay + 2V'{a) 



^1 - b+{a)/a - 2V{a) ^\ - b -{a) j a - 2V {a) j 
[(6+(a)/a)' + 2V'{a)f [(&-(a)/a)' + 2V\a)f 



4 \ [1 - 6+(a)/a - 2F(a)]3/2 [1 - 6_(a)/a - 2V{a)fl'^ 



(6+(a)/a)" + 2l^"(a) ^ (6_(a)/a)" + 2r'(a) 



Vl-&+(a)/a-2F(a) ^^1 - 6_(a)/a - 2y(a) j ' 



*V(«)\/i- 



fe+(a) 



-2y(a) + $'_(a)A/l- 



6_(a) 



- 2y(a) 



(63) 



(64) 



(65) 
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for which an easy computation yields: 

[47rS(a)a]' = + {$'|(a) ^1 - - 2l^(a) + Vl - h-{a)/a - 2y(a)} 

^f--(^) (^.(a)/ay + 2na) , ... (6-(a)/ay + 2^ (a) ] 

2\ + Vl -''+(«)/« -2^(a) " \/l-h-{a)/a-2V{a)\ 

and 

[47rS(a)a]" = {^'{'(a)^! - &+(a)/a - 2V(a) + ^'l{a)^l - h-{a)/a - 2y(a)} 

_ L„^^^ (6+(a)/a)^ + 2r(a) ^ ^„ (6_(a)/a)' + 2F'(a) | 



^l-6+(a)/a-2t/(a) " ^- («)/«- 2T/(a) J 

^ [(6+(a)/a)' + 2r(a)]^ ^ [(6„ (a)/a)' + 2^ (a)]^ 



[l-6+(a)/a-2y(a)]3/2 ' [I - b-{a) / a - 2V{a)]^'^ 



2l +^ ' ^l-h+{a)/a-2V{a) Vl - ^- («)/«- 2V'(a) J 



We shall evaluate these quantities at the assumed stable solution a^. 



H. The master equations 

In view of the above, to have a stable static solution at ao we must have: 



TOs(ao) 




6+(ao) 



6_(ao) 



ao 



'1 



ao 



(68) 



while 



and 



m'/(ao) > 



ms(ao) 
2ao 



l-y+(ao) ^ l-&L(ao) 1 
yi - 6+(ao)/ao a/1 - 6_(ao)/ao J 



(69) 



1 f [6+(ao) - ao6+(ao)]2 [6_(oo) - ao6'_(ao)]^ 



4ag 



[l-6+(ao)/ao]3/2 [1 - 6_(ao)/ao]3/2 



+ 



6;(ao) 



2 \ i/l - 6+(ao)/ao ^ ^Jl - b-{ao)/ao j" 



(70) 



This last formula in particular translates the stability condition y'(ao) > into a rather explicit and not too 
complicated inequality on m"(ao), one that can in particular cases be explicitly checked with a minimum of effort. 
For practical calculations it is more useful to work with ms{a)/a in which case 



[ms(a)/a]'| 



1 

+- 



(6+(a)/a)' 



+ 



(6_(a)/a)' 



(71) 



and 



1 



\n2 



[rUsW/al ^ +4 Ul-6+(a)/a]3/2 + [l-6_(a)/a]3/2 



1 

+ - 



(6+(a)/a)" ^ {b-{a)/a)" 



2 [ VI - b+{a)/a ^l-b-{a)/a) 



(72) 
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In the absence of external forces this inequahty (or the equivalent one for m"(ao) above) is the only stability constraint 
one requires. However, once one has external forces ^ 0), there is additional information: 



[AnE{a)a]'l^ = + {^'[{aWl - b+{a)/a + ^'^(a) ^1 - h.{a)/a]\ 



-y/1 - h+{a)/a 



Vl-6-(a)/a 



and (provided $±(ao) > 0) 

[47r5(a)an„„ < {^{0)^1 - 6+(a)/a + ^-(a)^! - &-(«)/«} | 



(&+(«)/«)' 

\/l - 5+(a)/a 



(6_(a)/a)' 



Vl-&-(a)/a, 



-l|$l(a)^iM^Ml^ + $' (a) 
4 r+^"^[l-6+(a)/a]3/2 + 

(6+(a)/a)" 



•\/l - b+{a)/a 



<i>'_(a) 



[(&-(«)/«)? 
[l-6_(a)/a]3/2 

(6_(a)/a)" 



^l-b-{a)/a 



If $j_(ao) < 0, we simply have 

[A-KE{a)a]"l^ > {$'|'(a)Vl - 6+(a)/a + $-(a) ^1 - &-(«)/«} | 



)/a]3/2 
(6+(a)/a)" , 



Vl-6_(a)/a 



^Jl - b+{a)/a 



<&'_(«) 



[(6-(«)/«)f 
[l-6_(a)/a]3/2 

(&-(«)/«)" 
Vl-&-(a)/a 



(73) 



(74) 



(75) 



Note that these last three equations are entirely vacuous in the absence of external forces, which is why they have 
not appeared in the literature until now. 



III. APPLICATIONS 



In discussing specific examples one now "merely" needs to apply the general formalism described above. Several 

examples are particularly important, some to emphasize the features specific to possible asymmetry between the two 
universes used in traversable wormhole construction, some to emphasize the importance of NEC nonviolation in the 
bulk, and some to assess the simplifications due to symmetry between the two asymptotic regions. 



A. Borderline NEC non-violation in the bulk: $± = 



On extremely general grounds we know the NEC must be violated somewhere in the spacetime of a traversable 
wormhole. By explicit computation we have seen that NEC violation in the bulk region external to the wormhole throat 
is equivalent to ^'±{r) < 0. Thus to minimize NEC violations for thin-shell wormholes we should demand ^±{r) > 
in the bulk. The boundary of this region corresponds to $± = (constant), which by a simple linear change of the time 

coordinates can be recast in the form $± = 0. Thus the situation $± = is particularly interesting not just because 
of mathematical simplicity, but physically interesting because it corresponds to the physical constraint that the bulk 
regions on either side of the wormhole throat be on the verge of violating the NEC. This $-|- ^ condition is satisfied, 
for instance, when the bulk spacctimes are Schwarzschild, Reissner-Nordstrom, de Sitter. Kottlcr (Schwarzschild- 
de Sitter), or Reissner-Nordstrom de Sitter, though the bulk spacetimes can be more general than any of these. Key 
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features of the $± = traversable wormholes are that in the bulk 

1 



p{r) = -prir) 



87rr2 



b', Mr) 



167rr 



b". 



while on the throat 



1 



1 - 



1 

Sira 



+ ■ 



(76) 

(77) 

(78) 



with the external forces vanishing (S = 0). Stability of an assumed static solution at ao then devolves into a single 
inequality (70) or the equivalent (72) being imposed upon m"(ao). 



B. Mirror symmetry: b± = b, $± = $ 

Another situation in which significant simplification arises is when the two bulk regions are identical, so that b± = b, 
and $± = Key features of these symmetric traversable wormholes are that on the throat 



(7 = — - 



27ra 



6(a) 



r = 



Ana 



aa - 



b{a)+ab' {a) 
2a 



with the external forces being rather simply determined by 

H = -$'(a) a. 



a 



(79) 
(80) 

(81) 



C. Thin-shell Schwarzschild wormhole: &± = 2M± and $± = 

Consider the situation where both bulk spacetimes on either side of the wormhole throat are portions of 
Schwarzschild spacetime. Then the metric functions of Eq. (1) arc given by 



&±=2M±, ^'±=0, 
respectively. For this case, inequality (70) yields the following inequality 



aoTO"(ao) > F(M±,ao) 



(M+/ao)^ 



(M_/ao)2 



(1 - 2M+/ao)'/' (1 - 2M_/ao)'/' 



(82) 



(83) 



The dimensionless function F{M±, qq) is depicted as the grey surface in Fig. 4, and the stability regions are situated 
above this surface. Wc have considered the definition x = 2Af+/ao for convenience, so as to bring infinite in to 
a finite region of the plot. That is, oq oo is represented as a; — > 0; and ao = 2Af+ is equivalent x = 1. Thus, 
the parameter x is restricted to the range < x < 1. We also define the parameter y ~ M^/M+, and from the 
denominator of inequality (83), wc verify that the parameter y lies within the range < y < 1/x. 

Prom a qualitative analysis of Fig. 4, we note that large stability regions exist for small values of x and of y. The 
stability regions decrease for large values of y, i.e., for M_ ^ A'/_|-, and for large values of x, i.e., for regions close 
to the event horizon. This analysis generalizes previous work, where the case of M_ = M+ was studied [28]. The 
stability regions decrease in the vicinity of the event horizon, a; ^ 1, and only exist for low values of y. 
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FIG. 4: Stability analysis for thin-shcU Schwarzschild traversable wormholcs. The stability region is that above the grey surface 
depicted in the plot. The grey surface is given by the dimensionless quantity F{Ad±,ao), defined by the right-hand-side of 
inequality (83). We have considered the range < a; = M+/ao < 1 and < y = M-/M+ < 1/x, respectively. Note that a 
large stability region exists for low values of a; = 2M+/ao and of j/ = M-/M+. For regions close to the event horizon, a; — >^ 1, 
the stability region decreases in size and only exists for low values of y. See the text for details. 



D. Thin-shell Reissner-Nordstrom wormholes 



The Reissner-Nordstrom spacetime is the unique spherically symmetric solution of the vacuum Einstein-Maxwell 
coupled equations. Its metric is given by 



^ 2M Q'^\ ,2 



(84) 



where M is the mass and is the sum of the squares of the electric (Qe) and magnetic (Qm) charges. In a 
local orthonormal frame, the nonzero components of the electromagnetic field tensor are F*"" = E'^ = Q-e/t^ and 
pecp — j^r _ Qy^jj-"^^ If |Q| < M, an event horizon is present, at a location given by 



n = M + - Q2. 



(85) 



If IQI > M we have a naked singularity. In the above, we have dropped the ± subscript as to not overload the 
notation. 

As in the Schwarzschild solution, one may construct a thin-shell wormhole solution, using the cut-and-paste proce- 
dure. For this case, inequality (70) yields the following dimensionless quantity 



aom;'(ao) > G(M±,0±,ao) = 



{Ml-Q\)/al 



(1 - 2M+/ao + Ql/alf'^ (l - 2M_/ao + Q''_/al) 



+ 



{Ml - Ql_)/al 



3/2 



(86) 
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FIG. 5: Stability analysis for thin-shcU Rcissncr-Nordstrom traversable worrriholcs. The stability region is that above the grey 
surface. The grey surface is given by the dirriensionless quantity G(A4±, Q±, ao), defined in Eq. (86). The values of \Q\/M+ = 1 
and \Q\/M+ = y/3/2 are depicted in the left and right plots, respectively. Note that for decreasing values of \Q\/M+, the 
stability regions decrease substantially for low values of M-/M+ and for high values of 2M+/oo. See the text for more details. 



The surface G{M±, Q±, ao) is shown in Fig. (5), and the stabiUty regions, ao 'm'siO'o), are depicted above this surface. 

As in the previous example, consider the following definitions: (i) x = 2M+/ao, in order to bring in infinity, i.e., 
ao oo is represented as x — >■ 0; and (ii) y = Using these definitions, and considering for simplicity that 

Q+ = Q- = Q, we have the following ranges 

< X < , < ?/ < ±^ . (87) 



^1 - Q^IMl 



The first range can be deduced from the condition that rb < ao < +oo; the second range can easily be deduced from 
the denominators of inequality (86). 

To now analyse the stabihty regions, we simply consider specific values for |Q|/M+. The values of |(5|/M+ = 1 and 
|(5|/M+ = •\/3/2 are depicted in the left and right plots, respectively, of Fig. 5. Note that for decreasing values of 
|Q|/M_|_, the stability regions decrease substantially for low values of M_/M+ and for high values of M+/ao. More 
specifically, if one were to construct thin-shell Reissner-Nordstrom wormholes with the junction interface close to the 
event horizon, one would need high values for the charge in order to have stable solutions. 



E. Thin-shell variant of the Ellis wormhole: 6± = It±/r 

It is perhaps instructive to consider an explicit case which violates some of the energy conditions in the bulk. 
Consider, for instance, the case given by the following shape functions: h± = R^/r, (so that b±(r)/r = R^/r'^), which 
were used in the Ellis wormhole [36]. (Note that Ellis' terminology is slightly unusual, and we have rephrased his 
work in more usual terminology.) 



1. Zero momentum flux: <1>± = 



As an initial step, consider the absence of external forces, so that 5 = 0. For this case one verifies that the NEC is 
borderline satisfied in the bulk. However, the WEC is violated as one necessarily has negative energy densities, which 
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FIG. 6: The plot depicts the stability region for the thin-shell variant of the Ellis wormhole, for which b± — R±/r. The grey 
surface is given by the dimensionless quantity H{R±, ao), given by the right-hand-side of inequality (89). The stability region, 
represented by the function oom"(oo), is that above the grey surfaces depicted in the plot. We have considered the range 
< a; = R+/ao < 1 and < y = R-/R+ < 1/x. Note that large stability regions exist for small values of x and of y. The 
stability regions decrease for large values of y, i.e., for R- ^ R+, and for large values of x. See the text for details. 



is transparent from the following stress-energy profile 



pir) = -pr{r) = -pt{r) = ■ (88) 



Stability regions are dictated by inequality (70), which yields the following dimensionless quantity 



aomg[ao) > H{H±,aQ) — + ^ — ^7372' ( 



The function H{R±,ao) is depicted as the grey surface in Fig. 6 and the stability regions lie above this surface. We 
consider the definition x = R+/ao for convenience, so as to bring in infinity, i.e., ^ oo is represented as a; ^ 0; 

and flg = Rj^ is equivalent to x = 1, so that the range for x is given by < x < 1. Wc also consider the definition 
y = so that from the denominator of H{R±, ao), given by Eq. (89), we verify that the range of y is provided 

by < y < 1/x. The stability analysis is similar to the Schwarzschild thin-shell wormhole case, in that large stability 
regions exist for small values of x and of y. The stability regions decrease for large values of y, i.e., for R- » 
and for large values of x. 



18 



2. Non-zero external forces: $± = —R±/r 

In order to generalize the previous specific example, we now consider a case with external forces. Specifically, let 
us consider the following functions 

$± = -— . (90) 
r 

These functions imply that $±(ao) > 0, so that in addition to the stabihty condition given by inequahty (70), 
one needs to take into account the stabihty condition dictated by (74). The latter inequality yields the following 
dimensionless quantity 

s,. > (6J^^■/»o - 19fiya§ + 12fl^/a§) , (6fl-/»o - IQfii/gg + 12fli/<.§) 

„.14.M„)l = (,.nx,.lf' ^ {^-Rlhlf ■ ' ' 

Note that as before, inequality (70) yields the dimensionless quantity given by inequality (89), which is depicted as 
the grey surface in Fig 6. The stability regions are depicted above this surface. We consider once again the definition 
X = R+/ao, so that the range for x is given by < a; < 1; and the definition y = so that as before the 

parameter y lies within the range; < y < l/x. 

Now, in addition to the imposition of inequality (89), depicted above the surface in Fig. 6, the stability regions 
are also restricted by the condition (91), depicted in the left plot of Fig. 7. In the latter, the stability regions are 
given below the grey surface. Collecting the results outlined above, note that the stability regions are given from the 
regions above the surface in Fig. 6, and below the surface provided by the left plot of Fig. 7. These stability regions 
are depicted in the right plot of Fig. 7 in between the grey surfaces, represented by the functions H{R±,ao) and 
H2{R±,ao), respectively. Note the absence of stability regions for high values of x and y. 




FIG. 7: The plot depicts the stability regions for the thin-shell variant of the Ellis worniholc, with b± = R±/r, but now in 
the presence of non-zero external forces. The latter non-zero momentum flux arises from the functions given by 3>± = —R±/r. 
The grey surface depicted in the left plot, is given by i?2(-R±,ao), i.e., the right-hand-side of inequality (91), and the stabihty 
regions are given below this surface. We have considered, as before, the range < a; = R+/ao < 1 and < y = R-/R+ < l/x. 
Now, collecting the results, note that the stability regions are given from the regions above the surface in Fig. 6, and below the 
surface provided by the left plot. These stability regions arc depicted in the right plot in between the grey surfaces, represented 
by the functions H{R±,ao) and H2{R±, ao), respectively. Note the absence of stability regions for high values of x and y. See 
the text for more details. 
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F. New toy model: b± = y/rR± 

1. Zero momentum flux: $± = 

Another new and interesting toy model is given by the following 

b± = ^/rR± and $± = 0, so that b±{r)/r = \jR±/r. (92) 
The bulk stress-energy profile is given by 

p(r) = -p,(r) = 4p,(r) = . (93) 

Note that the NEC and WEC are satisfied throughout the bulk. For this case, inequality (70) yields the following 
dimensionless quantity 



16(l-y^^)'^' 16(l-/Rl7^)^ 



aom"(ao) > /(i?±, ao) = — 7^ + — '^/^, ■ (94) 



The stability regions are presented in Fig. 8. As in the previous examples, we consider the definition x = R+/ao 
for convenience, so that the range for x is given by < x < 1. We also consider the definition y = R^/R+, so that 
the range of y is provided byO<y<l/x. Note that, as in the previous example of the thin-shell variant of the Ellis 
wormhole, large stability regions exist for small values of x and of y. The stability regions decrease for large values of 
y, i.e., for ::g> R+, and for large values of x. 




FIG. 8: The stability regions for the toy-model thin-shell wormhole with b± — ^rR±, <I>± = are given above the grey 
surfaces depicted in the plot. The grey surface is given by the dimensionless quantity I{R±,ao), given by the right-hand-side 
of inequality (94). We have considered the range < a; = R+/ao < 1 and < y = R-/R+ < 1/x. Large stability regions exist 
for small values of x and of y. Note that the stability regions decrease for large values of y, i.e., for R- ^ R+, and for large 
values of x. 
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2. Non-zero external forces: $± = —R±/r 

As in the previous application, it is useful to generalize the case of zero external forces. In order to have the presence 
of a non-zero momentum flux term, consider once again the following functions 

$±=-— , (95) 
r 

which imply that <I>±(ao) > 0. Thus, in addition to the stability condition given by inequality (70), one needs to 
take into account the stability condition dictated by (74). From the latter inequality, one deduces the following 
dimensionless quantity 



3u -(^V'^Tiu ^ {96R+/ao - 214(i?+/ao)3/2 + 117R%/al) 
a^[Ana^{a)\ >l2{R±,ao) = 

16 (l - /RTk) 



, (96i?-/ao - 214(fl_/flo)3/2 + n7R^_/al) 

+ . s 3/2 • 



16 [l-ylll/^y 



The stability regions arc dictated by inequalities (94) and (96). These arc depicted above the grey surfaces in Fig. 
8 and below the surface of the left plot in Fig. 9. The final stability regions are depicted in the right plot of Fig. 9 in 
between the grey surfaces, represented by the functions H{R±,ao) and H2{R±,ao), respectively. As in the thin-shell 
variant of the Ellis wormhole, note the absence of stability regions for high values of x and y. 




FIG. 9: The plot depicts the stability regions for the toy model thin-shell, with b± = \/R± r, but now in the presence of non-zero 
external forces. The latter non-zero momentum flux arises from the functions given by $± = —R±/r. The grey surface depicted 
in the left plot, is given by I2{R±, do), i.e., the right-hand-side of inequality (96), and the stability regions arc given below this 
surface. We have considered, as before, the range < x — R+/ao < 1 and < y — R-/R+ < 1/x. Now, collecting the results, 
note that the stability regions are given from the regions above the surface in Fig. 8, and below the surface provided by the left 
plot. These stability regions are depicted in the right plot in between the grey surfaces, represented by the functions I{R±,ao) 
and l2{R±,ao), respectively. Note the absence of stability regions for high values of x and y. See the text for more details. 
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G. Thin-shell chctrged dilatonic wormhole 

Consider a combined gravitational-electromagnetism-dilaton system [37-39] , described by the following Lagrangian 

L = V^{-i?/87r + 2(V?A)2 + e2*FV47r} . (97) 

In Schwarzschild coordinates the solution corresponding to an electric monopole is given by the following line element 
[37-39], 

ds' = -(l ^^=1 dt'+(l ^— 1 dr' + {de' + sin^ ^ V) , (98) 

where we have dropped the subscripts ± for notational convenience. The non-zero component of the electromagnetic 
tensor is given by F^f = Q/r'^, and the dilaton field is given by e^^ = 1 — Q"^ /M{j3 + + /32). The parameter (i is 
defined by /3 = QV2M. 

In terms of the formalism developed in this paper, the metric functions are provided by 



b{r) 




2M 



(99) 



$(r) =-l\n[l + '—\. (100) 



An event horizon exists at = 2M^Jl — fi/M. Note that $'(r) is given by 

which is positive throughout the spacetime. Thus, the stability regions are restricted by the inequalities (70) and 
(74). 

We consider for simplicity /3+ = The expressions for inequalities (70) and (74) are extremely lengthy, so we 
will not write them down explicitly. We define the following parameters x = 2M+/a and y = The junction 

interface lies within the range < a < oo, so that the range of the parameter x is given by 

< a; < ^ . (102) 

In the following analysis, we only consider the regions within the range < y < 1, as the stability regions lie within 
this range, as will be shown below. 

Consider as a first example the case for [3 ~ 1/2, so that the stability regions, governed by the inequalities (70) 
and (74), are depicted in Fig. 10. The left plot describes the stability regions above the grey surface, given by (70), 
while the right plot describes the stability regions below the grey surface, i.e., inequality (74). Note that the stability 
condition dictated by the inequality (70) shows that the stability regions decrease significantly for increasing values 
of a;, within the considered range of y, i.e., < y < 1. In counterpart, inequality (74) shows that the stability regions 
decrease significantly for decreasing values of y. These latter stability regions further decrease for increasing values of 
X. As a second case consider the value /? = 1/8, depicted in Fig. 11. We verify that the qualitative results are similar 
to the specific case of /3 = 1/2, considered above. 

Collecting the above results, the final stability regions, depicted in between the surfaces of different shades of grey, 
are given in Fig. 12. The left and right plot are given for the values of /3 = 1/2 and /3 = 1/8, respectively. Note that 
the stability regions decrease for decreasing values of (3. More specifically, for decreasing values of P, stability regions 
exist for practically very low values of y, i.e., M_ < Af+. In this context, it is interesting to note that in the vicinity 
of the event horizon, the stability regions increase for increasing values of j3, provided one has low values of y. 



IV. DISCUSSION AND CONCLUSION 



In this work, we have developed an extremely general, flexible and robust framework, leading to the linearized 
stability analysis of spherically symmetric thin shells. The analysis is well-adapted to general spherically symmetric 
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FIG. 10: Thin shell dilaton wormhole for /3 — 1/2. The left plot describes the stability regions above the grey surface, given 
by (70), while the right plot describes the stability regions below the grey surface, i.e., inequality (74). From the left plot, 
one verifies shows that the stability regions decrease significantly for increasing values of x, within the considered range of 
< y < 1. In counterpart, the right plot shows that the stability regions decrease significantly for decreasing values of y. These 
latter stability regions further decrease for increasing values of x. See the text for more details. 




FIG. 11: Thin shell dilaton wormhole for /3 = 1/8. The left plot describes the stability regions above the grey surface, given 
by (72), while the right plot describes the stability regions below the grey surface, i.e., inequality (74). The left plot shows 
qualitatively that the stability regions decrease significantly for increasing values of x, in the relevant range of < y < 1. From 
the right plot, one verifies that that the stability regions decrease significantly for decreasing values of y. See the text for more 
details. 



thin shell traversable wormholes and, in this context, the construction confines the exotic material to the thin-shell. 
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FIG. 12: The plots depict the final stabihty regions for the thin shell dilaton wormhole. The left plot describes the stability 
regions given by /3 = 1/2, while the right plot describes stability regions given hy P — 1/8. The final stability regions are 
depicted in between the surfaces of different shades of grey. Note that the stability regions decrease for decreasing values of /3. 
More specifically, for decreasing values of /3, stability regions exist for practically very low values of y, i.e., M- < M+. In this 
context, it is interesting to note that in the vicinity of the event horizon, the stability regions increase for increasing values of 
/?, provided one has low values of y. See the text for more details. 



The latter, while constrained by spherical symmetry is allowed to move freely within the bulk spacetimes, which 
permits a fully dynamic analysis. To this effect, we have considered in great detail the presence of a flux term, which 
has been widely ignored in the literature. This flux term corresponds to the net discontinuity in the conservation law 
of the surface stresses of the bulk momentum flux, and is physically interpreted as the work done by external forces 
on the thin shell. 

Relative to the linearized stability analysis, we have reversed the logic flow typically considered in the literature, 
and introduced a novel approach to the analysis. We recall that the standard procedure extensively used in the 
literature is to define a parametrization of the stability of equilibrium, so as not to specify an equation of state on 
the boundary surface [28-30]. More specifically, the parameter r]{a) = dV /da is usually defined, and the standard 
physical interpretation of 77 is that of the speed of sound. In this work, rather than adopt the latter approach, we 
considered that the stability of the wormhole is fundamentally linked to the behaviour of the surface mass rag (a) of 
the thin shell of exotic matter, residing on the wormhole throat, via a pair of stability inequalities. More specifically, 
we have considered the surface mass as a function of the potential. This novel procedure implicitly makes demands 
on the equation of state of the matter residing on the transition layer, and demonstrates in full generality that the 
stability of thin shell wormholes is equivalent to choosing suitable properties for the material residing on the thin 
shell. 

We have applied the latter stability formalism to a number of specific examples of particular importance: some 
presented to emphasize the features specific to possible asymmetry between the two universes used in traversable 
wormhole construction; some to emphasize the importance of NEC non-violation in the bulk; and some to assess the 
simplifications due to symmetry between the two asymptotic regions. In particular, we have considered the case of 
borderline NEC non-violation in the bulk. This is motivated by the knowledge that, on extremely general grounds, 
the NEC must be violated somewhere in the spacetime of a generic traversable wormhole, and if this were to happen 
in the bulk region, this would be equivalent to imposing '^'j^{r) < in the bulk. Thus to minimize NEC violations 
for thin-shell wormholes, we have considered the specific case of <i?± = 0, which is particularly interesting for its 
mathematical simplicity, and for its physical interest as it corresponds to the constraint that the bulk regions on 
either side of the wormhole throat be on the verge of violating the NEC. We have also considered the simplification 
when the two bulk regions are identical, and analysed the stability regions of asymmetric thin-shell Schwarzschild 
wormholes and thin-shell Reissner-Nordstrom wormholes in great detail. It was instructive to consider explicit cases 
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which violate some of the energy conditions in the bulk. For instance, we considered thin-shell variants of the Ellis 
wormhole and two new toy models, and explored the linearized stability analysis in the presence of zero momentum 
flux and non-zero external forces. Finally, we analyzed thin-shell dilatonic wormholes, where the exterior spacetime 
solutions corresponded to a combined gravitational-electromagnetic-dilaton system. 

In conclusion, by considering the matching of two generic static spherically symmetric spacetimes using the cut-and- 
paste procedure, we have analyzed the stability of spherically symmetric dynamic thin-shell traversable wormholes 
— stability to linearized spherically symmetric perturbations around static solutions. The analysis provides a general 
and unified framework for simultaneously addressing a large number of wormhole models scattered throughout the 
literature. As such we hope it will serve to bring some cohesion and focus to what is otherwise a rather disorganized 
and disparate collection of results. A key feature of the current analysis is that we have been able to include "external 
forces" in the form of non-zero values for the metric functions $±(r). (This feature is absent in much of the extant 
literature.) Another key aspect of the current analysis is the focus on ms{a), the "mass" of the thin shell of exotic 
matter residing on the wormhole throat — and the realization that stability of the wormhole is fundamentally linked 
to the behaviour of this exotic matter via a pair of simple and relatively tractable inequalities. 
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